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SUMMARY  PAGE 


PROBLEM 


To  improve  a  method  for  approximating  functions  of  wavelengths 
which  is  applicable  to  computer  calculations  related  to  vision. 

FINDINGS 

Shifting  the  origin  of  the  wavelength  scale  to  the  center  of  the 
region  of  interest  resulted  in  a  significant  reduction  of  roundoff  errors. 


APPLICATION 

Use  of  the  method  presented  enables  researchers  to  use  convenient 
approximation  methods  for  calculations  related  to  vision,  without  restriction 
on  the  computer  or  calculator  employed.  A  previously  published  method 
would  fail  on  some  machines  due  to  roundoff  error. 
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When  it  is  desired  to  fit  a  polynomial  to  a  function  of  wavelength,  roundoff  errors  may  be  minimized  by  shifting  the 
origin  of  the  independent  variable  (wavelength  or  frequency)  to  the  center  of  the  region  of  interest. 


As  mathematical  modeling  of  visual  functions  becomes 
more  successful,  it  is  naturally  desirable  that  independent 
workers  should  be  able  to  reproduce  a  given  model  easily  and 
accurately.  Boynton  and  Wisowaty1  have  taken  an  important 
step  in  this  direction  by  calculating  formulas  that  approximate 
a  certain  set  of  receptor-sensitivity  functions  put  forth  by 
Smith  and  Pokorny.2  In  this  Letter,  we  present  some  simple 
manipulations  of  the  Boynton  and  Wisowaty  results  that  will 
make  them  easier  to  use.  The  effect  of  the  transformation 
that  we  make  is  to  reduce  roundoff  errors,  but  the  goal  is  not 
to  make  the  answers  more  accurate.  Rather,  the  use  of  for¬ 
mulas  that  are  more  immune  to  roundoff  trouble  removes  the 
requirement  for  a  computer  that  carries  10  significant  digits 
and  the  need  for  coefficients  carried  out  to  10  digits. 

The  step  that  we  take  is  to  shift  the  origin  of  the  indepen¬ 
dent  variable  to  the  center  of  the  region  of  interest.  This 
traditional  device  of  applied  mathematics  has  general  appli¬ 
cation  when  a  polynomial  is  to  be  fitted  to  a  function  of 
wavelength.  If  the  device  were  applied  at  the  beginning  of 
a  regression  calculation,  it  would  reduce  or  eliminate  roundoff 
error  at  all  stages.  Finally,  in  the  same  spirit  of  making  cal¬ 
culations  more  workable  and  routine,  we  make  explicit  the 
relation  of  the  present  receptor  sensitivities  to  Judd’s  x',  y', 
and  z'  and  to  the  receptor  sensitivities  incorporated  in  the 
model  of  Guth  et  al .3 

In  discussing  the  use  of  their  formulas,  Boynton  and  Wi¬ 
sowaty  observe  that  “. . .  for  some  machines  rounding  errors 
may  cause  problems  with  the  fourth-order  polynomials  (it  is 
essential  to  retain  all  five  decimal  places  for  the  coefficients 
of  Table  I).”  4  By  “five  decimal  places,”  Boynton  and  Wi¬ 
sowaty  mean  some  8-10  significant  digits  on  most  of  the 
coefficients.  Roundoff  errors  occur  because  the  formulas  as 
given  depend  on  the  subtraction  of  nearly  equal  numbers.  For 
instance,  one  step  in  the  evaluation  of  the  green  receptor 
function  involves  the  polynomial  evaluation 

log  a  =  a4X4  +  a3X3  +  a^X2  +  a{\  +  a<j.  (1) 

Using  coefficients  oo  through  a4,  as  provided  by  Boynton  and 
Wisowaty,  and  X  =  0.5  fim  leads  to  the  sum 

log  a  =  -1747.361113  +6533.32095  -  9149.3564 
+  5695.5237  - 1332.9406 
=  -0.813463.  (2) 

To  shift  the  origin  of  wavelength,5  we  define  the  new  vari- 


As  a  practical  matter,  we  put  the  new  origin  in  the  middle  of 
the  visible  spectrum: 

X0  -  0.55  gm.  (4) 

Equation  (1)  may  be  replaced  by  a  polynomial  in  l: 

log  a  =  a4/4  +  a3l3  +  «2l2  +  otil  +  «o-  (5) 

In  general,  to  transform  the  coordinates  of  a  polynomial  of 
degree  N,  the  new  coefficients  ctk  can  be  calculated  from  the 
old  ak  by 

°tk-T.  ai+k  — rrr-  Mi,  (6) 

J—o  j!«! 

where  k  =  0. .  JV.  With  X0  =  0.55  fim  and  X  -  0.5  fim,  the 
evaluation  of  Eq.  (5)  involves  the  sum 

log  a  =  -0.1747  +  1.1551  -  2.7525  +  2.0055  -  1.0470 
=  -0.8136.  (7) 

The  answer  is  essentially  the  same,  but  roundoff  error  is  less 
likely  to  occur  since  we  are  no  longer  adding  and  subtracting 
quantities  104  times  as  large  as  the  answer.  Any  time  one 
wishes  to  fit  a  polynomial  to  data  that  have  a  limited  domain 
of  the  independent  variable  not  centered  on  the  origin,  a 
similar  reduction  of  roundoff  problems  can  be  expected  from 
a  transformation  such  as  Eq.  (3),  with  X0  in  the  center  of  the 
domain. 

In  the  case  of  the  red-  and  blue-sensitivity  functions, 
Boynton  and  Wisowaty  found  it  necessary  to  fit  polynomials 
in  frequency  to  functions  expressed  as  quantum  sensitivity. 
In  this  case,  the  origin  of  frequency  must  be  shifted.  We  de¬ 
fine 

f  =  v~v  0.  (8) 

As  a  practical  matter,  we  choose 

i’o  =  5.9  (9) 

in  units  of  1014  Hz — the  convenient  unit  selected  by  Boynton 
and  Wisowaty. 

The  working  equations  now  become  Eqs.  (3),  (4),  (8),  and 
(9)  above,  and  Eqs.  (10)— (22)  below.  Green-catching  sensi¬ 
tivity  is  given  by 

G  =  (an  +  bn)l/n,  (10) 

log  a  =  a4l4  +  a3(3  +  a2l2  +  ail  +  a0,  (11) 

log  b  =  /?4/4  +  fizl3  +  /?2/2  +  (12) 
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Table  1.  Values  of  a  for  Eq.  (10)  for  G,  and  for  Eq. 
(14)  for  R  and  B,  and  the  coefficients  for  Eqs.  (11), 
(12),  and  (15)— (17) 


Exponents 

R 

G 

B 

ni  -  2.25 

n  =  1.45 

ni  =  8.00 

n2  ~  8.50 

- 

n2  =  5.00 

Coefficients 

R 

G 

B 

a3  =  -0.02467 

a 4  =  -27960. 

a3  =  0.6877 

a2  =  -0.9920 

a3  =  -9241. 

<*2  =  -4.020 

Oil  =  1,432 

* — l 
O 
r— < 
r~i 

1 

II 

cq 

S 

«i  =  5.691 

ao  =  —2.087 

cn  =  -40.11 

cto  -2.053 

- 

a0  =  -1.047 

- 

=  0.5818 

04  =  -1706. 

03  =  -0.3913 

02  =  -1.042 

03  =  1190. 

02  =  -1.305 

01  =  -1.403 

02  =  -319.7 

0,  =  2.431 

00  =  -0.5700 

01  =  3.979 

0o  =  -0.7289 

- 

0o  =  -0.4519 

- 

73  =  0.6602 

- 

73  =  -11.05 

72  =  -0.6721 

- 

72  =  -33.31 

7i  =  -1.883 

- 

7i  =  -29.86 

7o  =  -0.9875 

- 

7o  =  -11.16 

Table  2.  Values  of 

n  and  Coefficients  for  Eqs.  (19)- 

(22)  To  Be  Used  in  Calculations  of  Eye  Lens  and 

Macular  Pigment  Density 

Exponents 

ni  =  15.00 
n2  =  7.50 
ri3  =  4.00 

Coefficients 


a3  -  -10130. 

73  =  -928.0 

a2  -  -3596. 

72  =  -485.1 

«!  =  -432.8 

7i  =  -56.18 

cro  =  -16.91 

7o  =  -1.295 

0s  =  -6684. 

63=  14.93 

02  =  -2243. 

d2  =  2.715 

0,  =  -244.5 

5i  =  -1.119 

0o  =  -8.016 

S0  =  0.0716 

Red-catching  sensitivity  R  and  blue-catching  sensitivity  B  are 
given  by 

R  or  B  =  Q (A/555),  (13) 

where  A  is  in  nanometers, 

Q  =  [(an>  +  bni)nRn'  +  cn2]C/"2);  (14) 

and 

log  a  -  a3f3  +  a2f2  +  a  if  +  a0,  (15) 

log  b  =  fo/3  +  /?2/2  +  Pif  +  do>  (16) 

log  c  =  73/3  +  72 f2  +  7i /  +  7o-  (17) 

Table  1  gives  the  required  constants. 


For  the  density  of  eye  lens  plus  macular  pigment, 


D  =  {[(oni  +  6ni)n2/ni  +  cn2]n3/n2  +  dni\lln*,  (18) 
where  the  equations  for  a,  b,  c,  and  d  are 

log  a  =  asl3  +  a2l2  +  ot\l  +  «o,  (19) 

log  b  =  fol3  +  P 2l2  +  Pi l  +  Po,  (20) 

log  c  =  y3l3  +  72(2  +  7i?  +  70,  (21) 

log  d  —  §3 13  +  d2l2  +  4-  §o-  (22) 

Table  2  gives  the  required  constants. 


The  polynomial  coefficients  have  been  tabulated  only  to 
four  significant  digits.  This  degree  of  precision  was  chosen 
in  a  systematic  way.  We  wrote  a  computer  program  to  tab¬ 
ulate  calculated  R,  G,  B,  and  D  with  the  polynomial  coeffi¬ 
cients  rounded  to  2,  3,  4, . . .  up  to  8  significant  digits.  Six¬ 
digit  coefficients  give  nearly  perfect  agreement  with  the 
Boynton  and  Wisowaty  calculated  values;  more  digits  give 
perfect  agreement.  Four-significant-digit  precision  in  the 
coefficients  gives  discrepancies  from  those  earlier  calculated 
values  of  a  few  units  in  the  fourth  decimal  place — a  precision 
of  about  1  in  104.  This  should  be  tolerable,  since  Boynton  and 
Wisowaty  find  their  fits  to  the  Smith  and  Pokorny  functions 
to  have  precision  of  the  order  of  2  in  102. 

SOME  RELATED  DATA 

Taking  the  inverse  of  the  transformation  matrix  derived  by 
Smith  and  Pokorny2  gives  the  expression  to  recover  Judd’s 
modified  x,y,  z6: 


(r) 

- 

f2.9448 

-3.5010 

0.2118\ 

r 

= 

1.0000 

1.0000 

0.0000 

G 

V'i 

0.0000 

0.0000 

1.0000 j 

This  provides  smooth  versions  of  the  Judd  functions  for  ap¬ 
plications  where  this  might  be  convenient. 

The  equations  above  give  (R,  G,  B)  with  maxima  of  (0.6442, 
0.3912, 1.6431)  at  wavelengths  of  (562,  542,  440)  nm.  Guth 
et  al.  utilize  essentially  the  same  (R,  G,B)  as  a  starting  point 
for  their  vector  model  for  normal  and  dichromatic  color  vi¬ 
sion,3  but  with  the  functions  adjusted  to  have  maximum  val¬ 
ues  of  1.0.  To  rescale  R,  G,  and  B  to  have  maxima  of  unity, 
multiply  the  expressions  given  above  [Eqs.  (10)  and  (13)]  by 
1.5523,  2.556,  and  0.6086,  respectively. 
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